Electron-phonon coupling and electron self-energy in electron-doped graphene: 
calculation of angular resolved photoemission spectra. 
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We obtain analytical expressions for the electron self-energy and the electron-phonon coupling in 
electron-doped graphene using electron-phonon matrix elements extracted from density functional 
theory simulations. From the electron self-energies we calculate angle resolved photoemission spec- 
tra. We demonstrate that the measured kink at « —0.2 eV from the Fermi level is actually composed 
of two features, one at « —0.195 eV due to the twofold degenerate E2g mode, and a second one at 
« —0.16 eV due to the mode. The electron-phonon coupling extracted from the kink observed 
in ARPES experiments is roughly a factor of 5.5 larger than the calculated one. This disagreement 
can only be partially reconciled by the inclusion of resolution effects. Indeed we show that a finite 
resolution increases the apparent electron-phonon coupling by underestimating the renormalization 
of the electron velocity at energies larger than the kinks positions. The discrepancy between theory 
and experiments is thus reduced to a factor of « 2.2. From the linewidth of the calculated ARPES 
spectra we obtain the electron relaxation time. A comparison with available experimental data 
in graphene shows that the electron relaxation time detected in ARPES is almost two orders of 
magnitudes smaller than what measured by other experimental techniques. 

PACS numbers: 74.70.Ad, 74.25. Kc, 74.25. Jb, 71.15.Mb 



Despite the fact that the band structure of graphene 
is calculated in many solid-state textbooks P, Q , its ex- 
perimental verification has been provided only recently 
by Angular Resolver Photo-Emission (ARPES) measure- 
ments on a graphene monolayer deposed on a SiC sub- 
strate The experiments show that the peculiar 
features of the electronic structure predicted theoretically 
are qualitatively confirmed by experiments: the carbon 
TT- bands (i) cross at the K-point in the Brillouin-zonc 
(Dirac point) and (ii) depart linearly with a slope Vf from 
the Dirac point, (iii) the Fermi velocity extracted from 
experiments d, Q is slightly larger (10-20%) than that 
calculated theoretically using density functional theory 
(DFT). 

Beside this encouraging agreement between theory and 
experience, recent angular photoemission (ARPES) ex- 
periments 0, 0] performed on graphene revealed remark- 
able surprises. Two kinks are seen in the ARPES dis- 
persion: the first one is at energies of 0.2 eV below the 
Fermi level (e/) and its energy position respect to e/ is 
unchanged as a function of the doping level while the 
second one is closer to the Dirac point and its energy po- 
sition respect to e/ decreases rapidly as the doping level 
is increased (see Fig. 2 in Ref. @). The first kink has 
been attributed to a phonon feature [6| , while the second 
kink has been interpreted as due to a plasmonj^, Q. In 
what follows we focus on the first kink. 

The ARPES momentum distribution curves (MDCs) 
associated to the -0.2 eV kink display a puzzling behav- 
ior as a function of doping. Indeed it is observed that the 
magnitude of the jump associated to the MDC-linewidth 
in the — 0.5eV < e — e/ < eV energy window decreases 
as a function of doping (see Fig. 3 in Ref. [6|, where 
from top to bottom the jump increases). This is surpris- 



ing since if this jump is associated to the electron-phonon 
interaction then it should refiect the imaginary part of 
the electron self-energy due to the electron-phonon inter- 
action. Since the magnitude of such interaction is usually 
proportional to the density of states at the Fermi level, 
the jump should increase as the doping level is increased. 
Thus the opposite behavior should be expected. This 
contradiction can be solved by noting that at low doping 
the tail of the second peak (attributed to a plasmon), is 
fairly close in energy and could affect the low energy part 
of the momentum distribution curve. At larger dopings 
the plasmon peak has no effect, the electron phonon 
coupling does increase as a function of doping. Thus we 
focus in the doping region identified by e/ > 0.3 eV (the 
energy-zero being at the Dirac point). 

In this work we calculate the electron-phonon coupling 
parameter and the electron-phonon coupling contribu- 
tion to the electron self-energy in doped graphene. In 
particular, we give an explicit demonstration of eq. 1 in 
ref. [l3|. From the electron self-energy we obtain the 
spectral-weight function and the ARPES spectra. Fi- 
nally we compare the calculated spectra with available 
experimental data, discussing in the details the impor- 
tant finite-resolution effects. 

The paper is structured as follows. In section |T| we 
obtain an analytical expression for the electron-phonon 
coupling in doped graphene. The electron self-energy is 
calculated in sec. |lll The ARPES spectra are calcu- 
lated from the electron self-energy in sec. IIIIl including 
finite resolution effects. Finally in sec. IIVI we compare 
the electron relaxation time measured by different exper- 
imental techniques for both electron-doped graphene and 
graphite. Sec. |V]is devoted to conclusions. 
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I. ELECTRON-PHONON COUPLING IN 
DOPED GRAPHENE 

In units of 27r/a with a — 2.46A, the 2D volume of the 
graphene Brillouin-zone (BZ) is = The graphene 

TT* bands are hnear with slope (3 — hvf = 5.52eVA 
(within DFT) close to the Dirac points K = (1/3, 1/3, 0) 
and K = 2K. The density of states per spin at a general 
energy e above or below the Dirac point, but still in the 
region where the tt* bands can be considered linear, can 
thus be written as 



27rV3|6 



47r|e| 



(1) 



Note that in this work energies are always measured re- 
spect to the Dirac point. 

The electron-phonon coupling for a mode v at momen- 
tum q due to the tt* bands in graphene is given by: 



A, 



huJciuNa{tf) JbZ ^ 

S{ek - e/)(5(ek+q - £/) 



dk ^ 

Q I5k7r*k+q7r* 



(2) 



where is the phonon frequency of the mode at mo- 
mentum q and ffk^.k+qTr* electron-phonon matrix 
element for the 7r*-band ek and for the phonon mode v. 

To illustrate how the integral is evaluated we introduce 
the following two regions of space, namely the sets: 

TK{e) = {k|/3|k-K| < e + ry} (3) 
.FK'(e) = {k\f3\k-K'\<e + rj} (4) 

In these definitions, is a small positive quantity. For 
e = ef, since we assume that the Fermi level is not too 
far from the Dirac point so that the tt* bands are lin- 
ear, |k — K| or |k — K I is a small but finite vector and 
•^K (e/ ) n -^K' (^/ ) is empty. The boundary of each region 
of space at r/ = (circumference) is indicated as dJ-ji{ef) 
and dT^' (e/). 

In Eq. [21 the two (5— functions restrict the k inte- 
grations to the region of space satisfying the conditions 
Ek = £/ and Ek+q = e/- The set of k points such that 
fk = £/ is composed by the set dJ^ic{£f){JdT^'{ef). 
Thus in the integral in Eq. [2l two cases are given (label- 
ing k' = k + q): 

(i) k,k' ea.FK(e/) or k,k' G a.fK'(e/), 

(ii) k e dTi^{ef) , k e d^K.'i^s) O'' '^i*^^ versa 

In case (i) scattering occurs at q = F + q, with small q, 
and TT* bands can only couple to the twofold degenerate 
E2g phonon mode. In case (ii) scattering occurs at q = 
K-fkoratq=K +k, with small k , and tt* bands can 
only couple to the Aj^ phonon mode [11[. The electron- 
phonon matrix elements involved in the two scattering 



process have been fitted to ah initio data in Ref. |ll| 
and are: 



^K+k7r*,K+k+q7r* ' 
'K+k7r%K'+k+q7r' ' 



(5^)[l±cos( 



'k.q 



k,k+q 



a)l5) 



(.4) 



1 + cos( 



k,k+q^ 



(6) 



In Eq. [5] the ± sign refers to the LO/TO E2g modes 
respectively, {gl) = 0.0405eV2 and (g^.) = 0.0994 eV^, 
and 0u,v is the minimal angle between the two vectors 
u,v. 

For case (i) one has 



qB2 



2 X 2 X 2{gl)F f (fk 

X^(eK+k - e/)'^(eK+k+q - e/) 



In 



where 



TT{e) = {-k\l3k<e + v}. 



(7) 



(8) 



The 8 prefactor is the results of having 2 E2g modes and 
of having an identical integral over the second Fermi sur- 
face sheet at K'. The integral /q is the so-called nesting 
factor, defined as: 

f (Pk 

h= -F^'^l^K+k - e/)^(eK+k+q - e/) (9) 

The electron-phonon coupling due to E2g modes is 
given by 



Ar(e/)-A^(6/) / ^Aq^,^ 



2N{ef){g. 



hwTE2 



where we have used that 



^r(2e^) 



(10) 



(11) 



and we have replaced the E2g phonon frequency with its 
value at T. 

Similarly, case (ii) leads to 



A 



2 X 2{g^^)F 



(fk 



X'^(eK'+k-e/)^(«K+k+q-e/) 

■^K + q 



[1 - C0S(% ii+q)] X 



and the additional factor of 2 is a result of having scatter- 
ing from jrK(e/) to.?^K'(^/) and vice versa. The quantity. 



(12) 



<fk 



X'^(eK'+k-^/)'5(^K+k+q-^/) (13) 



3 



EC k, co)= 




FIG. 1: Lowest order contribution to the electron self-energy 
due to the electron-phonon interaction. The dotted (continu- 
ous) line represents the phonon (electron) self-energy. 



and its integral over the momentum q are evaluated in 
sec. IVII A[ so that the contribution of the mode to 
the electron-phonon coupling is: 



AK(e/) 



(14) 



KA, 



where we have approximated <^K+qyi' ~ ^k.a' ■ 
The total electron phonon coupling is thus: 



X{ef)^N,{ef) 
which is eq. 1 in ref. 
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FIG. 2: (Color online) Real ( E ) and imaginary (E ) parts of 
the electron self-energy in graphene (continuous line). Dashed 
lines refer to self-energy parts otained using a constant density 
of states. The Fermi level is e/ = 0.4 eV 



II. ELECTRON SELF-ENERGY AND ANGLE 
RESOLVED PHOTOEMISSION 

The lowest contribution to the retarded electron self- 
energy due to coupling of tt* electrons to a phonon mode 
V is illustrated in Fig. [TJ At zero temperature direct 
calculation [13] of the diagram gives: 



S.(k,e) 



q, ,2 ^ 

~Pp|yk7r'.k-Fq7r' I ^ 
BZ " 



E 

={-1,1} 

6(ae/ - aek+q) 



_e + i5 - (ck+q) + aftwq 



(16) 



where Q{x) is the Heaviside function. The imaginary 
part of eq. [H] is 



BZ 



I5k7r*,k+q7r* I 



X e(Q:e/ - aek+q)^(e - ek+q + aftwq 
= — TT 9(Q!e/ — ae — ?iajqiy) 

a={-l,l} 



(17) 



In angular resolved photoemission (ARPES) experi- 
ments the graphene is electron-doped, so the Fermi level 
is larger than the Dirac point but it is still in the region 
where ek can be considered linear. For a given mode u 
and a given value of a, the 5-function in eq. 1171 restricts 



the BZ integration to two regions, close to K and to K . 
The restriction to these regions of k-space and the fact 
that we are interested in the region of ener_gy-momentum 
close to the Dirac point, namely k = K -|- k with k small, 
restricts furthermore the integration region. Indeed it 
implies that for small q: (i) q = q and (ii) q = K -|- q. 
Case (i) represents scattering to phonons close to the F 
point while (ii) to phonons close to the K point. So the 
situation is similar to the previous electron-phonon cal- 
culation. 

The total self-energy, E — X]iy={£;2 A'}^iy^ "^^^ 
the two phonon modes at T and to the Ai phonon 
mode at K is obtained substituting Eqs. [5] and [5] in Eq. 
1171 assuming a constant phonon dispersion around F and 
K and performing the integration over the BZ, as 

"={-14} 

x9(ae/ — ae — /iwrBsg) + ^KA'^-^K(e — "^kAj) ^ 

X 6(ae/ — at — hjj-^^i 



(18) 



where Ar(e — hjjTE2g) ^^nd AK(e — ^"^ka' ) ^'^^ defined in 
Eq. [To] and in Eq. [Ml respectively. From Eq. [18] we 
note that for small k the imaginary part of the phonon 
self-energy is momentum-independent, so in what follows 
we drop the k-label. Using numerical values of hLUrE2g — 
0.195eV and ^^ka' ~ 0.16eV, the S (e) is illustrated in 
Fig. [1 (black hues)'. 

The imaginary part in eq. [TH]has to be compared with 
the square well model which is obtained from Eq. [TS] 
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assuming a constant density of states. This is the com- 
monly used approximation to interpret ARPES spectra 
[H, m. The square well model is illustrated in fig. [2] 
(red-dashed) . In graphene this approximation is in prin- 
ciple not allowed due to the behavior of the density of 
states proportional to |e| (see eq. [T]) . The difference 
between the two models becomes relevant for energies 
smaller than or closer to the Dirac point. 

The real part of the electron self-energy can be ob- 
tained using the Kramers-Kronig relations, namely: 



E(e) 



E (6) 



de 



(19) 



If the self-energy in Eq. [18] is used then E (e) diverges 
at large |e| due to the |e| dependence of the density of 
states. However this divergence is unphysical since in real 
graphene the linearity of the tt* bands and the consequent 
behavior of the density of states is only up to energies 
of « 1.5 eV from the Dirac point. Thus for large |e| 
the imaginary part of the phonon self-energy should be 
regularized. We adopted the following regularization: 



■'Reg 




if |e| < en 
if EM < |e| 



(20) 



With this assumption, 
eq. [19] leads to: 



the calculation of the integral in 



E(e) 



+ 



~N4ef)(gl) {{e 
log|(eA/ + e)(e/ 

(e - ncjrB^ J log 



- hL0TE2 



[-2 log I £ 
-)|] + 



e/ -I- hujTE^g - e 



{< 



loe 



2 

(em + e){ef -hiu 



e + fi^KA, ) -2 log 



+ 



(e- SwKA , )log 



•' 1 



(gf )(e - hwTE2„) + -Tr{<^ - fii^KA, 



log 



EM + e 



EM 



where we allowed for a small constant imaginary part rj 
to eliminate numerical instabilities. Typically -q = 10""* 
eV 



III. RESULTS 

In this section we consider a Fermi level of e/ = 0.4 eV 
measured from the Dirac point. Neglecting resolution ef- 
fects, the spectral function (Eq. [22]) is shown in figure 
[3] (left) . The scans presented are fixed photoelectron en- 
ergy scans (the energy value is given on the left of Fig. 
[3] (left) ) while the photoelectron momentum is varied 
(MDC scans). In the MDC scans two kinks are present 
(compare scans g and h for the first kink and scans e and 
d for the second), both in the spectral- weight maximum- 
position and linewidth. 

The behavior of the MDC-maximum-position as a 
function of energy and momentum is illustrated in fig. [3] 
(right). The lower energy kink corresponds to the twofold 
degenerate E2g mode, while the higher energy one corre- 
sponds to the mode. 

The behavior of the MDC-linewidth, half-width half- 
maximum (HWMH), as a function of energy is shown 
in fig. U] Notice that, in absence of resolution effects, 
the linewidth is equal to — /3E (e). The behavior of the 
linewidth is compared with that of square well model 
typically used to interpret ARPES data. The differences 
are negligible at energies larger than -0.2 eV, but they 
are significative at lower energies and in particular at the 
Dirac point. 

To test the robustness of the phonon features against 
experimental resolution we introduce the following con- 
voluted spectral weight: 



de / d^k 



Ae,p(k,e) = /(e) 



A(k ,e')G^,(6' -e)G^Jk - k) (23) 

where G,,^ (cc) is a Gaussian having full- width rjx and cen- 
tered in a; = 0. The Fermi distribution is indicated with 

1 



/(e) = 



exp[(e-e/)/fcBr] + l 



(24) 



In practical calculations we choose em = 1.5 eV. 

The real part of the electron self-energy is illustrated 
in Fig. [2] and is compared with the real part obtained 
from the Kramers-Kronig transformation of E {lo) with 
a constant density of states. 

In ARPES experiment the spectral weight is measured, 
namely 



This form of the experimental resolution assumes to have 
decoupled momentum and energy resolutions. 

(21) 

We chose % — 25 meV and 77k = 0.1 A , as in recent 
ARPES experiments [1, Q. The maximum position in 
Aoxp(k, e) is plotted in Fig. [5] 

As it can be seen the experimental resolution substan- 
tially smears out the two kinks so that a very weak kink 
is visible at -0.2 eV while the second one is almost invisi- 
ble. The kink is substantially smaller than what detected 
in experiments [a| . The mass enhancement parameter is 
defined as [13, [iJ: 



A{k,e) 



-2[E"(e)+7y] 



[e-ek-E'(e)]^-f [E"(e)+^?]■ 



(22) 



5E (e) 



de 



(25) 
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FIG. 3: MDC scans for different energies (left) and ARPES maximum position as a function of energy and momentum (right). 
The Fermi level is e/ = 0.4 eV, the experimental resolution is not included in the calculation. A 1 meV broadening is used for 
illustration purposes. 



Linearizing E (e) 
(for ?7 = 0): 

A(k,e): 

where 



-Ae, the spectral-weight becomes 



-2E"(e)Z2 



[e-Zek]' + [S"(e)Z]' 



Z = 



1 + A 



(26) 



(27) 



is the quasiparticle weight. From Eq. [26] one sees that 
the quasiparticle state has quasiparticle energy Zck and 
linewidth S (e).Z/2. In graphene the bands are linear, 
with Efc = /3A:, so that the maximum position in the spec- 
tral weight at energies higher than the kink is given by 
the relation 



1 + A 



(28) 



Assuming linear rcnormalizcd bands, ek — /3phfc, for en- 
ergies larger than the kink then the following expression 
for A is obtained: 



A = 



- 1 



(29) 



Typically is obtained from a linear fit to maximum 
position in MDC curves at energies below the kink (in our 
case /? = 5.52eVA within DFT), while /3ph is obtained 
from a linear fit at energies higher than the kink but 



enough below e/ so that the effects of the Fermi function 
in eq. [21] are absent. 

When a finite resolution is used, the result is substan- 
tially affected. Indeed, we find that a linear fit in the en- 
ergy window from — 0.195 eV < e — < — 0.04eV leads 
to values of the electron-phonon coupling which are a fac- 
tor 2.5 than Eq. [TS] as shown in Fig. [5] A similar fitting 
procedure is used in experiments (note that in experi- 
ments only the kink due to the Yi2g phonon is visible). 
A comparison with available experimental data (see fig. 
5 in Ref. 0) is shown in fig. [H] A clear disagreement 
between theory and experiments is still present, even if 
it is reduced by resolution effects. 

For this reason, we believe, the determination of A us- 
ing Eq. [29] is affected by a large error. The error is 
due to the difhculties in the determination of /3ph gen- 
erated by the non-linearity in E (e) near the kink. The 
non-linearity, when convoluted with a finite resolutions, 
results in a quasi-linear behavior with an "apparent" en- 
hanced electron-phonon coupling. 



IV. ELECTRON RELAXATION TIMES IN 
GRAPHITE AND GRAPHENE 

A. Electron-doped graphene 

It is interesting to compare the ARPES-measured self- 
energy imaginary part with what detected by alternative 
experimental techniques. In electron-doped graphene the 
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FIG. 4: (Color online) MDC - linewidth (HWMH) for the 
two model self-energy. The curve in upper panel is obtained 
using a finite resolution, while that in the lower one is with- 
out resolution effects. In the absence of resolution effects the 
linewidth is equal to — /3E (e) (see fig. [TJ . The Fermi level is 
ej = 0.4 eV. 
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FIG. 5: Position of the maximum in MDC curves using a 
finite resolution. The Fermi level is e/ = 0.4 eV. 



is performed (T w 300 K ). As shown in appendix I VII Bl 
this leads to electron scattering time of the order of 

T = G.35ps (from mobility) (30) 

for £/ > 0.2eV. Since the Debye temperature of the op- 
tical phonon in graphene is much larger then 300 K, the 
scattering in the mobility measurements is mainly due to 
defects and acoustic phonons. 

ARPES measures the electron self-energy at photo- 
emitted electron-energies e. The imaginary part of the 
electron self-energy is related to the electron scattering- 
time by the relation T(e) — h/[2Yj (e)]. From the data in 
Refs. [1,01 for |e — e/| ~ ksT we obtain 

T « 3.5fs (from ARPES) (31) 

which is two order of magnitudes smaller than Eq. 1301 



electron relaxation time has been determined experimen- 
tally by conductivity /mobility data [Iql and by angular 
resolved photoemission measurements [q. The mobility 
measurement detect the scattering time of electrons at 
energies |e — e/| < fc^T, where ks is the Boltzmann con- 
stant and T is the temperature at which the experiment 



B. Graphite 

In graphite the electron scattering time has been 
measured by two different experimental techniques; (i) 
Femptoseconds time-resolved spectroscopy [l3l and (ii) 
ARPES. 
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From femptoseconds time- resolved spectroscopy [17[: 

T w 0.2ps for |e- e/l = 0.25 eV 
r w O.lps for |e- e/l = 0.50eV 
(from femtoseconds 

time — resolved spectroscopy) (32) 

Similar to what happens in graphene, ARPES mea- 
surements [3, H m, [3 lead to a relaxation time which 
are two order of magnitudes times smaller than what 
obtained from the femptoseconds photoemission spec- 
troscopy. For example, from the measured ARPES 
linewidth in Ref. [l^ (see Fig. 10 c), we obtain: 

T « 4.7fs for |e-e/| =0.25eV 
r « 3.4fs for |e- e/| = 0.50 eV 

(from ARPES) (33) 



and is of the same order of magnitude of the electron 
relaxation time obtained from conductivity and fem- 
toseconds time-resolved spectroscopy measurements in 
graphite. However this is in strong disagreement with 
ARPES measurements, being the ARPES relaxation 
times, both in graphene and graphite, almost two or- 
der of magnitudes smaller. This discrepancy essentially 
reflects the disagreement in the measured and calculated 
electron-phonon coupling. The aforementioned disagree- 
ment in the electron self-energies is even more surprising 
when considering that previous DFT calculations of the 
phonon self-energy in graphene, graphite and nanotubes 
were found to be in perfect agreement with experimental 
data for what concerns phonon dispersion and phonon 
lifetimes, [U |2^. Since the electron and phonon 
self-energies involve the same vertex, and thus the same 
matrix elements, a good agreement would be expected 
even for the electron self-energies too. 



V. CONCLUSIONS 



In this work we calculated the electron-phonon cou- 
pling parameter and the electron-phonon coupling con- 
tribution to the electron self-energy in doped graphene. 
From the electron self-energy we obtained the spectral- 
weight function and the ARPES spectra. 

The ARPES spectra as a function of momentum and 
energy displays two kinks. The kinks are at energies e — 
Cf « —0.195 eV and e — e/ ~ —0.16 eV, where e/ is 
the Fermi level. The two kinks are due to coupling to 
the twofold degenerate mode and to the mode 
respectively. The MDC-linewidth as a function of energy 
is discontinuous (jump) at E2g and A^ phonon energies. 

Comparing the calculated electron-phonon coupling 
with that extracted from ARPES experiments we found 
that, for large enough electron-doping, the latter is 
roughly a factor 5.5 larger than the former, as suggested 
in ref. [oj. We partially solved this contradiction by in- 
cluding finite resolution effects. Indeed, in experiments 
the electron-phonon coupling is determined from the ra- 
tio of the electron-velocities at higher and lower energies 
respect to the kink. The velocities are obtained from the 
slopes of the maximum position of the ARPES spectra 
as a function of energy and momentum. We find that 
the slope above the kink is substantially affected by the 
presence of a finite resolution and the extracted values 
of the electron-phonon coupling are « 2.5 larger than 
what obtained without any resolution effect. Thus when 
comparing calculated spectra with the inclusion of finite 
resolution effects to ARPES experiments M remark 
that the measured electron-phonon coupling is still a fac- 
tor of 2.2 larger than the calculated one Thus this 
work shows once more [20.] the importance of including 
resolution effects to correctly describe ARPES data. 

Finally, from the imaginary part of the electron self- 
energy we obtain the electron-relaxation time. The 
calculated electron relaxation time is in good agree- 
ment with mobility data on electron-doped graphene 
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VII. APPENDIX 

A. Evaluation of the integral Jx+q 

In Eq. m %,k+q = 20kq-'^' ^o that l-cos(% £+-) = 
2cos^(6'j^-) and the integral Jk+q (see Eq. [T^] ) can be 
evaluated as: 



xS{f3\k- 



dk ^~ 



2kcos^{e^^)d{(3k- (3kf) x 



2ir 



^Si^^kj + 2A:gcos(6'j^-) 
^ Y: £ d(cos(0g.))5(cos(0ej + _f ) X 



kf) 



2k f 



cos 



kfq\ sin(0e.) 



q/2kf 



r 



(34) 
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Furthermore, the integral of such a quantity is: 
9 r 



1 



2ir 2 

d6i / dy- ^ 



arcsin(2/) 



167rfc2 _ 16^1 _ Nl{tf) 



(35) 



B. From mobility to electron relaxation time 



The conductivity tensor is [2J| 

= 2e'E/ (02^-(kK(k)-v,(k) 



5/ 
(9e 



(36) 

«=e.(k) 



Eq. should then divided by two since we are interested 
in only one component of the conductivity tensor. More- 
over we only consider intraband transition since they are 
the only relevant at low temperature. The graphite bands 
are linear so that = ±/3fc = ±hvfk where k measure 
the distance from the K— point. At zero temperature 
only electrons on the Fermi surface contribute to the inte- 



gral thus, assuming a constant relaxation time, one gets: 

dkk '■2'^ 



(2^) 



d96{ef ~ Eki.) 



^2y-^£/ 

u J 



so that 



e^ef 



(37) 



(38) 



The conductivity can be written as a function of the mo- 
bility as (T = defj,, where S is the number of electrons 
participating in conduction per surface area. Using Eq. 
1 in Ref. 11 (5 = ej/iirf]^) and cr = ejefi/ {nf3'^). Thus 
the relaxation time becomes: 



(39) 



where £/ is expressed in eV and ^ in cm2/(Volt x sec) 
Expressing the mobility in m2/(Volt x sec) one gets 



r = 1.47 X 10 €fn = 1.47 x e//ips 



(40) 



Using the values of Ref. [l6| for mobility one gets at large 
doping (e/ > .18 eV) values of the order of 0.35ps. 
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